/// What do we owe soldiers? ///
/// Blagden, Krebs, and Ralston ///

* MS Table 2: Observational, Knoweldge of the Covenant
reg favor milcovkno  if check==1  & control ==1
estimates store m1, title(Model 1)
reg favor milcovkno milkno ideo miltherm hawk female age education income nonwhite selfserv household milkno if check==1  & control ==1
estimates store m2, title(Model 1)
logit binarysensitivity  milcovkno if check==1 & control==1
estimates store m3, title(Model 1)
logit binarysens milcovkno  ideo miltherm hawk female age education income nonwhite selfserv household milkno if check==1  & control ==1
estimates store m4, title(Model 1)
esttab m1 m2 m3  m4   using "UK Cov Table 2 Paper.rtf", se r2 pr2  replace starlevels(+ 0.10 * 0.05 ** 0.01 *** 0.0001) b(2) se(2)


* Figure 1
logit binarysens i.milcovkno  if check==1  
margins i.milcovkno
marginsplot, recast(bar)

* Figure 2
logit milcovkno  ideo miltherm hawk female age education income nonwhite selfserv household milkno if check==1  
estimates store m5, title(Model 1)
coefplot m5, scheme(s1mono) mlcolor(black) mfcolor(none) mcolor(black) ///
		xline(0, lpattern(dash) lcolor(red%50)) ///
		drop(_cons) mlabposition(6)  ///
		aspect(1) xlab(-2.5 "Less knowledge"  0  2.5 "More knowledge")

* MS Table 3: Experimental, Prime of the Covenant
reg favor milcovprime  if check==1  
estimates store m1, title(Model 1)
reg favor milcovprime  ideo miltherm hawk female age education income nonwhite selfserv household milkno if check==1  
estimates store m2, title(Model 1)
logit binarysensitivity milcovprime if check==1 
estimates store m3, title(Model 1)
logit binarysens milcovprime   ideo miltherm hawk female age education income nonwhite selfserv household milkno  if check==1  
estimates store m4, title(Model 1)
esttab m1 m2 m3  m4   using "UK Cov Table 3 Paper.rtf", se r2 pr2  replace starlevels(+ 0.10 * 0.05 ** 0.01 *** 0.0001) b(2) se(2)

* Figure 3
reg favor i.milcovprime##c.ideo if check==1
margins milcovprime, at(ideo=(0(0.2)1))

marginsplot, ///
recast(line) ///
recastci(rarea) ///
ciopts(color(%20)) ///
plot1opts(lcolor(blue)) ///
plot2opts(lcolor(red)) ///
legend(order(1 "Control" 2 "Covenant")) ///
title("Interaction Effect Plot with 95% CI") ///
xtitle("Ideology") ///
graphregion(color(white))
